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Abstract 

The possibility to verify the pseudo-Dirac nature of neutrinos is inves- 
tigated here via the detection of ultra high energy neutrinos from distant 
cosmological objects like GRBs. The very long baseline and the energy range 
from ~ TeV to ~ EeV for such neutrinos invokes the likelihood to probe very 
small pseude-Dirac splittings. The expected secondary muons from such neu- 
trinos that can be detected by a kilometer scale detector such as ICECUBE 
is calculated. The pseudo-Dirac nature, if exists, will show a considerable 
departure from flavour oscillation scenario in the total yield of the secondary 
muons induced by such neutrinos. 
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1 Introduction 



Evidence has been obtained from the satellite-borne observations, the ex- 
istence of the Gamma Ray Bursts (GRB) from extra galactic (or galactic) 
sources characterised by sudden intense flashes of 7-rays. The GRBs can pro- 
duce very high energy (> 10 16 eV) cosmic rays. Although the mechanism of 
such GRBs are not fully understood but various model calculations suggest 
that neutrinos of that energy range are also produced in GRBs and should 
be detected by very large detectors like a kilometer cube water cherenkov 
detector (e.g. ICECUBE at south pole). Because of their astronomically 
long baseline (^ Mega Parsec) these ultra high energy (UHE) neutrinos may 
open up a new window in very small mass square difference (Am 2 ) regime 
and may throw more insight to neutrino physics. A viable way to detect such 
UHE neutrinos is to look for upward going secondary muons produced by the 
charged current (CC) interaction of such UHE with the terrestrial rock. 
Already installed AMANDA in the south pole and the future ICECUBE [1] 
detector, are able to detect such secondary muons. Both are water Cerenkov 
detectors and use south pole ice as the detector material. 

The possibility that the UHE neutrinos from a distant GRB can probe 
very low Am 2 region much below the solar or atmospheric neutrino regime, 
can be utilised in investigating the proposed pseudo-Dirac nature of the neu- 
trinos [2, 3] through its detection in large (~ Km 2 ) detectors [4] like ICE- 
CUBE. In the pseudo-Dirac scenario, each active neutrino is split into two 
almost degenerate components of active and sterile part. This can be theo- 
retically realised from the generic mass matrix in [ul, {vr) c ] basis which can 
be written as (for one generation) 



where the tiny Majorana mass terms and itlr (m^,m^ << m^) are 
introduced to slightly lift the degeneracy of Dirac mass m D . Thus the Dirac 
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neutrino is split into a pair of almost degenerate Majorana neutrinos with 
nearly maximum mixing angle given by tan 28 = 2m£,/(mf i — mi,). In this 
scenario, therefore, each of the three neutrino species is split up into a pair 
of neutrinos with very tiny mass square difference Am 2 = 2mn(mL + ttir). 
For three generation scenario, therefore each of the three types of neutrinos 
has this pseudo-Dirac splitting of very small mass square difference. 

The purpose of this work is to demonstrate the possibility of probing the 
pseudo-Dirac nature by studying the neutrino oscillation effects through the 
detection yield of neutrino induced muons in a large (Km 2 ) neutrino detector 
such as ICECUBE. For this purpose, the expected muon signal induced by 
neutrinos from a GRB at Mpc distance is calculated by separately folding (i) 
the flavour oscillation effects and (ii) the oscillation effects in pseudo-Dirac 
scenario to the GRB neutrino flux. 

2 GRB flux, neutrino oscillations and num- 
ber of secondary muons 

The GRB neutrino flux is estimated considering the the relativistic fireball 
model [5]. In the relativistic expanding fireball model of GRB, protons (also 
electrons, positrons and photons) produced in the magnetic field of the rotat- 
ing accretion disc around a possible black hole are accelerated perpendicular 
to the accretion disc at almost the speed of light forming a jet which is re- 
ferred to as fireball. The burst is supposed to be the dissipation of kinetic 
energy of this relativistic expanding fireball. These very highly energetic 
protons in the jet then interact with the photons and produce pions (cosmo- 
logical beam dump) through the process of A resonance. (The pions are also 
produced through the pp process). These pions then decay to yield and 
v e in the approximate proportion of 2:1. 

The neutrino flux from a GRB depends on several GRB parameters. 
Firstly, the Lorentz boost factor T is required for the transformation from 
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the fireball blob to observer's frame of reference. As the shocked protons in 
the blob photoproduce pions, the photon break energy (as the photon spec- 
trum is considered broken) and the photon luminosity L 7 (generally ~ 10 53 
ergs/sec) are important parameters for determining the neutrino spectrum. 

With all these, the neutrino spectrum from a GRB can be parametrised 
as [6, 7] 

^ = Axnmi(U/£;)xi (1) 
In the above, E u is the neutrino energy N v is the number of neutrinos and 

r 2 

E* ~ 10 6 b 2 - 5 GeV 

-^7,MeV 

r 2 . 5 = r/io 2 - 5 

A = E GRB 

1 + \n(E vmax /E b u ) ' 1 ; 

where E umilx is the cut-off energy for the GRB neutrinos and ^grb is the total 
energy that a GRB emits. Now, the observed energy E° bs of a neutrino with 
the actual energy E v coming from a GRB at a redshift distance z is given 
by the relation E° hs = E v /(1 + z) and similarly, the maximum observable 
neutrino energy E°^ x is E°^ x = E umax /(1 + z). The comoving distance d 
of a GRB at redshift z is given by 

d(z) = - / . (3) 

where Qm is the matter density (both luminous and dark), Q\ is the dark 
energy density respectively in units of critical density of the universe and c,H 
are the velocity of light in vacuum and Hubble constant respectively. In the 
present calculation c = 3 x 10 5 Km/sec and H = 72 Km/sec/Mpc (1 Mpc 
= 3.086 x 10 19 Km). Therefore the neutrinos from a single GRB that can be 
observed on earth per unit energy per unit area of the earth is given by, 

dN° hs dN v 1 



dE° hs dE v 4nd 2 (z) 



(1 + *) (4) 
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Here, as is evident from the above, the total flux per unit area per unit energy 
from a single GRB is being considered rather than the flux per unit time. 

The production process of UHE neutrinos suggests that the neutrino 
flavours are produced in the ratio v e : : v T = 1 : 2 : 0. Considering 
the maximal mixing between and v T (as indicated by the atmospheric 
neutrino data) and the element XJ e j, of the mass-flavour mixing matrix to 
be zero, the flavour ratio on reaching the earth, for neutrino mass-flavour 
oscillation, becomes v e : : v T — 1 : 1 : 1 irrespective of the solar mixing 
angle. Needless to say, because of the astronomical baseline (L ~ Mpc), the 
acquired relative phases of the propagating neutrino mass eigenstates is aver- 
aged out (Am 2 L/E » 1) and the UHE neutrinos from a GRB reaching the 
earth are incoherent mixture of mass eigenstates. Therefore, the probability 
for measuring a particular flavour j3 by a terrestrial neutrino telescope, if 
only flavour oscillation is considered, is Pp = 1/3. 

On the other hand, in pseudo-Dirac scenario, we have each of the three 
mass eigenstates vi, u 2 and u 3 to be nearly degenerate pairs and thus one 
obtains a total of six mass eigenstates. Kobayashi and Lim [3] worked out the 
mixing in such scenario and calculated the oscillation probability. Follwing 
[3] and Beacom [4] the probability Pp to detect a flavour (3 by a neutrino 
telescope for pseudo-Dirac neutrinos is 



P/3 = Y, w *H\ U «3\ 2 \ U Pi 



2 Ami 

1 - sin 2 ±- 

\ 4E 



(5) 



where m,j(j = 1,3) denotes the mass eigenstates for the three types of neu- 
trinos, Am 2 is the mass square difference due to pseudo-Dirac splitting of 
the mass eigenstate isj. a, (3, .. denote the flavour index and U a j is the CKM 
matrix for three generation mass to flavour mixing. w a is the relative flux of 
each of the neutrino flavours (a) at the production point (J2a w a — !)■ 

The total number of secondary muons induced by GRB neutrinos at a 
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detector of unit area is given by (following [8, 9, 6]) 

r E obs A /Vobs 

s = JfJ" d K hs ^PsUK hs , o z )P,(K hs , b^) (6) 

In the above, P SUTV is the probability that a neutrino reaches the detector 
without being absorbed by the earth. This is a function of the neutrino- 
nucleon interaction length in the earth and the effective path length X(9 Z ) 
(gm cm -2 ) for incident neutrino zenith angle 9 Z [9 Z = for vertically down- 
ward entry with respect to the detector). The interaction length L int is given 
by i 

L ' mt = a tot (Ef s )N A ^ 

and 

P surv (P° bs , 9 Z ) = exp[-X(9 z )/L iQt ] = exp[-X(9 z )a tot N A }. (8) 

where N A (= 6.022 x 10 23 gm _1 ) is the Avogadro number and a tot (= a cc + 
a NC ) is the total cross section. The effective path length X(9 Z ) is calculated 

as 

X(6 Z ) = J p(r(9 z ,£)d£. (9) 

In Eq. (9), p(r(9 z ,£) is the matter density inside the earth at a distance 
r from the centre of the earth for neutrino path length £ entering into the 
earth with a zenith angle 9 Z . The quantity P M (P° bs , E thr ) in Eq. (6) is the 
probability that a secondary muon is produced by CC interaction of and 
reach the detector above the threshold energy E thr . This is then a function 
of v^N (N represents nucleon) - CC interaction cross section a cc and the 
range of the muon inside the rock. 

P^(P° bs , E thl ) = N A a cc (R(E: hs ; P thr )> (10) 

In the above {R(E° hs ; E thl )) is the average muon range given by 

1 r l-E thl /E v rlrT CC (F, ohs 11) 

<P(P° bs ; E thr )) = -^ jf dyR(E^(l - y), E thr ) ™ ^ ' V) (11) 
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where y = (E° bs — E^)/E^° s is the fraction of energy loss by a neutrino of 
energy E° hs in the charged current production of a secondary muon of energy 
Ep. Needless to say that a muon thus produced from a neutrino with energy 
E v can have the detectable energy range between E t h_? and E v . The range 
R(En, E thr ) for a muon of energy is given as 

m>Bthi) f E 'J^^lJ^±^) (12) 

JE thI {dE^/dX) (3 \a + f3E thT ) 

The average lepton energy loss with energy E^ per unit distance travelled is 
given by [8] 

(§) =—^ (13) 

The values of a and (3 used in the present calculations are 

a = {2.033 + 0.077 ln[£ M (GeV)]} x 10~ 3 GeVcm 2 gm _1 

(3 = {2.033 + 0.077 ln[E M (GeV)] } x lO^Wgrn- 1 (14) 

for E^ < 10 6 GeV [10] and 

a = 2.033 x lO^GeVcmVn- 1 

P = 3.9 x 10" 6 cm 2 gm- 1 (15) 

otherwise [11]. 



3 Calculations and results 

The GRB neutrino flux is calculated for a GRB with energy -Eqrb — 10 53 
ergs. The neutrino break energy E h v is calculated following Eq. (2) with the 
Lorentz factor T = 50.12 and corresponding photon break energy E^ MeV = 
0.794. These values are obtained from Guetta et al [12] from their fireball 
model framework calculations. These values are tabulated in Ref. [6]. The 
calculation is performed for several values of redshift (z). In this calculation 
all length units are taken in Km. 
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The probability Pp in Eq. (5) is computed with solar mixing angle #© = 
32.31°, the atmospheric mixing angle # atm = 45.0° and 1-3 mixing angle 
013 = 0. Also, the representative pseudo-Dirac splittings are considerd as 
Am 2 = 10~ 12 eV 2 for each of the three species. Note that, in this case, Pp 
for mass-flavour oscillation is 1/3. GRB neutrino flux with redshift z, per 
square kilometer on earth per unit energy (GeV) for a zenith angle Z , is 
obtained using Eq. (1 -4) and then multiplied by the probability Pp. 

The secondary muon yield at a kilometre scale detector such as ICECUBE 
is calculated using Eqs. (6 - 15). The earth matter density in Eq. (9) 
is taken from [9] following the Preliminary Earth Reference Model (PREM). 
The interaction cross-sections - both charged current and total - used in these 
equations are taken from the tabulated values (and the analytical form) given 
in Ref. [13]. In the present calculations E°^ x = 10 11 GeV and threshold 
energy E thr = 1 TeV are considered. 

Fig. 1 shows the total yield of secondary muons induced by UHE neu- 
trinos from GRBs of different red shift values ranging from 0.02 to 3.8, for 
three cases namely (i) no neutrino oscillation (ii) mass-flavour oscillation and 
(hi) the psedo-Dirac oscillation. The zenith angle is taken to be 6 Z = 180° 
(vertically upwards neutrino, i.e. muons coming from vertically below the 
detector). Fig. 2 is same as Fig. 1 but for 9 Z = 100.9°. The pseudo-Dirac 
case distinctly differs from flavour oscillation and no oscillation scenario. In 
the present case considered here, the secondary muon events for pseudo-Dirac 
oscillation is almost half the yield from expected events if the UHE neutrinos 
suffers only mass-flavour oscillation and same is around six times less in case 
there is no oscillation. In Fig. 3, the secondary muon events in each bin of 
E° bs are shown. The events below 10~ 15 are not shown in Fig. 3 for clarity. 
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4 conclusion 



In this report, the possibility of probing pseudo-Dirac neutrino through ultra 
high energy neutrinos from distant GRBs is considered. It is demonstrated, 
that such UHE neutrinos from GRB are indeed capable of producing oscil- 
lation effects clearly distinct from flavour oscillation or no oscillation of such 
neutrinos. This can be realised by looking at the yield of secondary muons, 
produced by these neutrinos, at a kilometer scale detector such as ICECUBE 
at south pole. Detecting UHE neutrinos is also a viable means to probe as 
also verify the pseudo-Dirac splitting if it really exists. 
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Figure Captions 



Fig. 1 The comparison of total secondary muon yields in cases of (i) pseudo- 

Dirac oscillation (ii) Mass-Flavour oscillation and (iii) No oscillation for ultra 

high energy neutrinos from GRBs at different zshift (z) distance. The zenith 

angle 9 Z is 180°. See text for other GRB parameters considered. 

Fig. 2 Same as Fig. 1 but for 6 Z = 100.9°. 

Fig. 3 The secondary muon events in each bin of E° hs . 
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